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ABSTRACT 

We provide complete and global solutions of viscous transonic flows around 
black holes. We show that for any degree of advection, there are two critical 
viscosity parameters a c i, C 2 such that for a < a ci the flow may pass through the 
inner sonic point only. For a ci < a < a c2 , the flow may have shock waves, and 
for a > a c2 , the flow may again pass through the inner sonic point, depending 
on flow parameters. No new topologies emerge other than what we found earlier 
while studying viscous isothermal transonic flows. These findings corroborate 
our unification scheme discussed in the context of the stable, transonic solutions 
around black holes and neutron stars. 
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1. Introduction 



The study of viscous, transonic flows was initiated by Paczyhski & Bisnovatyi-Kogan 
(1981) to incorporate proper treatment of accretion disks close to a black hole, where radial 
velocity is significant, and where the standard Shakura-Sunyaev (1973) disks are likely 
to be unstable. This model was followed up by Abramowicz et al. (1988) where local 
computation of ra(E) including advective flux showed stabilizing effect on the inner region 
of the Shakura-Sunyaev disks. First global solution of the viscous, transonic flows were 
obtained by (Chakrabarti, 1990ab; hereafter C90ab). For simplicity, these works assumed 
that the heating and cooling were adjusted in a manner that the temperature of the disk 
remained constant. In the language of Shakura-Sunyaev (1973) viscosity parameter a, we 
showed that if a < a cr} the incoming flow may either have a continuous solution passing 
through outer sonic point, or, it can have standing shock waves at x s3 or x s2 (following 
notations of C90ab). For a > a cr} a shock wave at x s2 persisted, but the flow now had two 
continuous solutions — one passed through the inner sonic point, and the other through 
the outer sonic point. Later analytical and numerical works (Chakrabarti & Molteni, 1993, 
1995; hereafter CM93 and CM95 respectively) showed that x s3 is stable, and that for 
a > a cr the continuous solution passing through the inner sonic point is chosen. We noted 
that a cr (~ 0.015 for the case considered) was a function of the model parameters, such 
as the disk temperature, sonic point location and angular momentum at the inner sonic 
point. Our result did not agree with the findings of Muchotrzeb (1983), who also studied 
isothermal flows and found that solutions should not exist for some a > a* ~ 0.02. 

Subsequently, extensive numerical simulations of quasi-spherical adiabatic accretion 
flows (Molteni, Lanzafame & Chakrabarti, 1994), showed that shocks form very close to 
the location where vertically averaged model of adiabatic flows predict them (Chakrabarti, 
1989; hereafter C89). The flow advected its entire energy to the black hole and the entropy 
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generated at the shock is also totally advected allowing the flow to pass through the inner 
sonic point. It was also found, exactly as predicted in C89, that flows with positive energy 
and higher entropy form strong, supersonic winds. In presence of viscosity also, very little 
energy radiates away (e.g., Fig. 8 of C90a). Having satisfied ourselves by the stability 
of these solutions (CM93, MLC94, CM95), we proposed an unified scheme of accretion 
disks (Chakrabarti, 1993, 1994, 1995; hereafter C939495, CM95) which combines the 
physics of formation of sub-Keplerian disks with and without shock waves depending on 
viscosity parameters and accretion rates. We always considered only the stable branch of 
the transonic flow and our solutions remained equally valid for black hole and neutron star 
accretions as long as different inner boundary conditions are employed. 

Our present interest to revisit the problem arose after we became aware of newly 
discovered advection dominated models of accretion flows around black holes (Narayan & 
Yi, 1994; hereafter NY94) and the 'unified' scheme of accretion disks (Chen et al., 1995; 
hereafter CALNY95) which use cooling mechanisms different from C90ab. These solutions 
supposedly discuss new branches which become advection dominated close to the black 
holes. These solutions did not find shock waves (which were found to be very robust in 
simulations of CM93, CM95, MLC94, Ryu et al, 1995). Since no self-similar model may 
have any preferred length scales including sonic points and shock waves, shocks are not 
expected in NY94 treatment. Similarly, local solutions with the assumption of Keplerian 
angular momentum distribution = flRepierian as in CALNY95 could not have shock waves 
either, as the flow does not encounter any centrifugal barrier that is needed to form a shock 
wave (unless the flow is pre-heated by external radiations as in spherical accretion flows of 
Chang k Ostriker, 1985). 

In this Letter, we present exact global solutions of the equations used by Narayan 
& Yi (1994), (but applicable to black hole accretions) and show that no new topology 
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of solutions emerge other than what we had already obtained in the context of viscous, 
transonic, isothermal flows (C90ab). We neither assume the flow to be self-similar as in 
NY94 (and therefore Mach number of our disk is not a constant), nor assume the flow to be 
Keplerian as in CALNY95. As a result, as in isothermal flows (C90ab), we obtain solutions 
which allowed the disk to have shock waves as well. We also obtain a critical a parameter 
with exactly the same meaning as before. We therefore suspect that the so-called advection 
dominated flows do not constitute any new solutions. Our suspicion becomes stronger when 
we find that the conclusions of NY94 work regarding the advection of energy and entropy, 
formation of winds with positive energy and similarity of quasi-spherical results with those 
of vertically averaged model are exactly same as what we already discovered in references 
cited above. 



We assume units of length, velocity and time to be x g = 2GMbh / c 2 , c and 2GMbh / c 3 
respectively. In order to obtain complete and global topologies of the 'advection dominated' 
flows, we employ here the same set of equations as in NY94. 

(a) The radial momentum equation: 



2. Model Equations 
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(b) The continuity equation: 





(c) The azimuthal momentum equation: 
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(d) The entropy equation: 
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Here I Key 1S the Keplerian angular momentum, W, S and W X( f, are the vertically 
integrated density p and viscous stress respectively, h(x) is the half-thickness of the 
disk at radial distance x obtained from vertical equilibrium assumption (C89), i? is the 
radial velocity, s is the entropy density of the flow, Q + and Q~ are the heat gained and 
lost by the flow. Unlike NY94, we use Paczyhski-Wiita (1980) potential to describe the 
black hole geometry. We also drop self-similarity assumption (unlike Chakrabarti, 1990c, 
where self-similarity in vertically averaged flows were used for the first time in studying 
nox-axisymmetric flows). We employ Shakura-Sunyaev (1973) viscosity prescription 
W x< f, = aW } where W is the integrated pressure P. As noted earlier (CM95), it is probably 
more appropriate to use W X( f, oc II while studying shock waves since the total pressure 
II = W + S^ 2 is continuous across the shock and such a W X( f, keeps the angular momentum 
across of the shock to be continuous as well. Except for Eqn. Id, other equations are the 
same as used in our previous studies (C89 and C90ab). In C89, Eqn. Id was replaced by 
the adiabatic equation of state P = Kp 1 with entropy constant K different in pre-shock 
and post-shock flows, and in C90a, Eqn. Id was replaced by the isothermal equation of 
state, W = K 2 Yj (K being the sound speed of the gas). In the present Letter, for simplicity, 
we assume /(x)=const as in NY94. We have verified that our conclusions do not change 
when more general cooling laws are used instead. Details will be presented else where. 



3. Results 
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We solve equations l(a-d) using sonic point analysis (Liang & Thompson, 1980; C89, 
C90ab). Figs, l(a-b) show variations of Mach number M = $/a } (a 2 = ^fP/p } 7 being 
the polytropic index chosen to be 4/3 here) with the radial distance as a is varied. The 
angular momentum at the inner edge of the disk and the 'cooling parameter' / are chosen 
to be l e dge = 1-65 (as in C89) and / = 0.5 for illustration. The Figures are drawn with 
sonic points at x c = 2.8 in la and x c = 34.58 in lb respectively, exactly as in Fig. 7 of C89. 
For a = 0, the solutions are represented by the solid curves (as in Fig. 7 of C89). In Fig. 
la, as a (marked on the curves) is increased, the closed topology opens up (long dashed 
curves) till a > a c2 ~ 0.0975 is reached when the flow completely opens up (short dashed 
curve) to connect with Keplerian disks at larger distance. The location at which this joining 
takes place depends on a > a c2 , provided other parameters (x C} ledge, f) are kept fixed. In 
Fig. lb, as a is increased, open solutions (long dashed curve which could join Keplerian 
disks at a large distance) become closed (short dashed curve) when a > a c2 ~ 0.01. Other 
branch of the solution is also shown by a short dashed curve. Note that at the sonic point, 
M/ 1, but a value appropriate for our one-and-a-half dimensional model (C89). An exact 
expression for M(x c ) is: M 2 (x c ) = [a 2 f + n + 1 + {a 2 f(a 2 f + 1) + (n + l) 2 } 1 / 2 ]/[ 7 (2n + 1)] 
where, n = (7 — For a = 0, this goes over to M(x c ) = 2n/(2n + 1) as in C89 and over 

to M(x c ) = 1 for 7 = 1 (isothermal flow) as in C90a. 

We note here that the above consideration of critical a c2 applies only if the flow 
parameters are already in the 'accretion shock' region (i.e., where the outer sonic point has 
less entropy that the inner sonic point, C89). If the flow parameters are in the 'wind region' 
(i.e., when the inner sonic point has less entropy that the outer sonic point, or if the outer 
sonic point does not exist at all), then the flow will continue to pass through the inner sonic 
point (and join with the Keplerian disk) untill anothe critical viscosity parameter a ci is 
reached. 
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A solution may include a shock wave provided Rankine-Hugoniot conditions (or, 
any other generalized shock conditions, see, Abramowicz & Chakrabarti, 1980; C90b) 
are satisfied somewhere between the inner and the outer sonic points and the entropy 
at the inner sonic point is higher compared to that at the outer sonic point. Maximum 
a for which shocks could form with inner sonic point at x c is a c2 . As noted earlier, 
a Xc depends strongly on other parameters. For example, for x c = 2.8, and l e d ge = 1-65, 
a cr = 0.01, 0.015, 0.075, 0.0975, 0.165 for / = 0, 0.1, 0.3, 0.5, and 1.0 respectively. 
Figs. 2(a-b) show an example of the solution of Eqns. l(a-d) which includes a shock 
wave. These solutions were completely missed by NY94, Abramowicz et al. (1988) and 
CALNY95. Here, we choose a = 0.05 and / = 0.5 for illustration, and employ W X( f, = —all 
prescription to guarantee continuity of angular momentum through the shock. The solid 
curves represent the branch passing through the outer sonic point located at x out = x c = 50 
and the long dashed curves represent the branch passing through the inner sonic point 
at Xi n = x c = 2.8695. Initially, the flow passes through x out and becomes supersonic. 
After a standing shock at x s3 ~ 13.9 (vertical short-dashed curve) it becomes subsonic 
and enters the black hole through X{ n along the long dashed curve. The flow chooses 
this solution and not that by solid curve for x < x s3 since the entropy of the flow is 
higher at the inner sonic point (C89, CM93, MLC94). In Fig. 2(a), we plot the shock 
invariant (C89) function C, angular momentum distribution l(x) (x2), the Mach number 
distribution M(x), the total pressure II (in arbitrary units) and the local energy density 
E(x) = l/2$ 2 + (7 - l) _1 a 2 + l/2l 2 /x 2 - l/2(x - I)" 1 . The solid and dashed curves 
intersect at x = x s3 consistent with a Rankine-Hugoniot condition assuming an infinitesimal 
shock thickness. In Fig. 2(b), we show how the temperature T = /im p a 2 /~fk (in units of 
10 11 K) and the radial velocity 1? (in units of c) vary in a complete solution. Here, ji = 0.5 
for pure hydrogen, m p and k are proton mass and Boltzman constant respectively. At the 
shock, temperature goes up and the velocity goes down, in the same way as in our earlier 
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studies (C89, C90ab). We note here that for 7 > 1.5 the vertically averaged flow does not 
have three sonic points and therefore shock solutions are possible only if the flow is already 
supersonic. 

Since out transonic solutions with / =const have one less free parameter (as the flow 
satisfies sonic point conditions), the M(E) relation is monotonic with positive slope and 
represent the stable solutions. Whether or not a transonic solution is possible at all depends 
upon the nature of d$ / dx\ Xc (x C} a, /, l e d ge )- Here, we have left out the topologies which 
produce winds (with and without shocks) but we have verified them to be similar to the 
isothermal solutions (C90ab). Details will be presented else where. 

4. Discussion and Conclusions 

In this Letter, we presented complete and global solutions of viscous, transonic 
equations which may or may not include shock waves, depending upon the accretion rate 
(determining the cooling rate), angular momentum at the inner edge as well as viscosity 
parameters. Though a new name ('advection dominated', NY94) has been recently coined 
to describe a few of these solutions, we find no new topologies other than those already 
present in the literature in the context of viscous, transonic, isothermal flows (C90ab). 
Unlike self-similar solutions of NY94, which allowed only constant Mach number flows, our 
solutions show very rich topological properties in the M — X plane and also allow formation 
of shock waves. 

Upon discovering no new topologies with different cooling laws, we are once again 
convinced of our unified scheme (C939495, CM95) for flows around a black hole. If the 
accretion rate and viscosity are such that a shock wave can form, the post-shock solution 
would behave as a thick accretion flows discussed earlier in the literature (e.g., Paczyhski 
& Wiita, 1980; Rees et al., 1982) only more consistent, as advection is included as well. If 
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the viscosity is high enough, the flow becomes so-called 'slim' disk or 'advection-dominated' 
disk, only more self-consistent, as we neither assume self-similar flows (NY94) nor locally 
Keplerian flows (CALNY95). Furthermore, our angular momentum distribution does not 
show inexplicable undulations seen in 'slim' disks (Abramowicz et al, 1988), but smoothly 
joins with Keplerian, even when shock waves are present. We show that for a given cooling 
parameter (e.g., /), a solution can match with Keplerian flows at larger distance, either 
when a < a ci or if a > a c2 when the flow had three sonic points. If a d a < a c2 , the 
solution includes a shock wave. We know of no other solutions in the Literature which are 
as exact and complete as we discussed in C90ab, and in this Letter. Just as C89 and C90ab 
solutions were found to be stable through time-dependent numerical simulations (CM93, 
MLC94, CM95), we believe that our present solutions containing shocks are also stable. 
Indeed, we have already verified that disks with a power law (in T) cooling (which includes 
bremsstrahlung) have stable shock waves (Molteni, Sponholz & Chakrabarti, 1995). These 
shocks show oscillations (with 5 — 10 percent of amplitude in emitted radiations) whenever 
the cooling time roughly agrees with the in fall time and have all the characteristics of 
Quasi-Periodic Oscillations. Similarly, diversion of flows through the outer and inner sonic 
points mediated by viscosity (CM95) may be responsible for dwarf novae type outbursts. 
In another development, Chakrabarti & Titarchuk (1995) show that Comptonization 
of the soft photons (coming from the pre-shock Keplerian disks) in the hot post-shock 
region produces hard X-rays and 7-rays compatible with the observations of galactic and 
extra-galactic black hole candidates. 

Although, we have used pseudo-Newtonian potential, the result may be exactly 
similar when full general relativity is used. For instance, inviscid flows in Kerr geometry 
(Chakrabarti 1990d) show no new topology other than what is in C89. Even fully general 
relativistic studies of magnatohydrodynamic flows (Takahashi et al, 1990, Englemaier, 1993) 
show no new topologies other than what was obtained with psudo-Newtonian geometry 
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(Chakrabarti, 1990e). We therefore believe that our description of black hole and neutron 
star accretion is complete. 
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Fig. 1. — (a-b): Variation of Mach number with logarithmic radial distance as viscosity 
parameter a is varied. The advection parameter / = 0.5 is chosen. Solid curves (a = 0) 
represent solution of C89. At the inner sonic point (a), for a < a cr ~ 0.0975, and at the 
outer sonic point (b), for a < a cr ~ 0.01 (solid and long dashed curves), the topologies are 
useful for shock formation. For a > a cr} (short dashed curve) the topologies are open and 
useful for complete global solutions in (a). In (b), the topology becomes closed and this 
branch is not available for a black hole accretion. 



Fig. 2. — (a-b): An example of a complete solution which includes a shock wave, a = 0.05 
and / = 0.5 are used. Solid curves show quantities as matter passes from subsonic to 
supersonic at x out = 50. Shock conditions are satisfied at x = 13.9 (vertical dash curve) and 
the flow subsequently enters the black hole along the long dashed curves becoming supersonic 
at Xi n = 2.8695. See text for details. 



